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Abstract. 

In an irrotational dust universe, the locally free gravitational field is covariantly 
described by the gravito-electric and gravito-magnetic tensors E a b and H a b- In 
Newtonian theory, H a j, — and E a b is the tidal tensor. Newtonian- like dust universes 
in general relativity (i.e. with H a b = 0, often called 'silent') have been shown to be 
inconsistent in general and unlikely to extend beyond the known spatially homogeneous 
or Szekeres examples. Furthermore, they are subject to a linearization instability. Here 
we show that 'anti-Newtonian' universes, i.e. with purely gravito-magnetic field, so 
that E a {, = / Hab, are also subject to severe integrability conditions. Thus these 
models are inconsistent in general. We show also that there are no anti-Newtonian 
spacetimes that are linearized perturbations of Robertson- Walker universes. The only 
E a b = / H ao solution known to us is not a dust solution, and we show that it is 
kinematically Godel-like but dynamically unphysical. 



PACS numbers: 0420, 9880, 9530 



1. Introduction 

Irrotational dust spacetimes are characterized by vanishing pressure (p = 0)| and 
vorticity (u a = 0), and positive energy density (p > 0). They are important arenas 
for studying both the late universe and gravitational collapse models ||. Apart 
from the spatially homogeneous and isotropic Friedmann-Lemaitre-Robertson- Walker 
(FLRW) case, which is characterized by vanishing shear (cr ab = 0), these spacetimes have 
non-zero shear and non-zero locally free gravitational field. This field is represented by 
the irreducible electric and magnetic parts E a b and H a b of the Weyl tensor (see £§] for a 
covariant analysis of local freedom in the gravitational field). 

X It is implicit that the anisotropic stress Tr a b and the energy flux q a also vanish. 
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The gravito-electric tensor E a b is the relativistic generalization of the tidal tensor 
in Newtonian theory, while the gravito-magnetic tensor H a b has no Newtonian analogue 
|5|, P|, and is associated with gravitational radiation. There are no sound waves 
in dust, and when H a b = there can be no gravitational radiation |7|, §]. Dust 
spacetimes with H a b = have therefore been called 'silent', since there are no 



propagating signals [H], fnj.§ Within the class of irrotational dust models, the 
silent (H ab = 0) universes may be called 'Newtonian-like' , since they have a clear 
Newtonian counterpart. One might therefore expect that there is a broad variety of 
such Newtonian-like spacetimes, representing the general relativistic generalization of 
simple Newtonian models. However, there are subtleties involved in the Newtonian limit 
of general relativity (see for example || 112, 13|, Q); an d it turns out that the relativistic 
Newtonian-like spacetimes form a very restricted class. 

Independent exact analyses of silent H a b = models, using different methods, 
are given in |15| and ||16|| . The independent approaches produce the same conclusion. 
The silent condition H a b = leads to an integrability condition whose repeated time 
differentiation forms in general a non-terminating chain that leads to inconsistencies. It 
is thus unlikely that consistent silent solutions exist beyond the known special cases of 
Szekeres spacetimes and spatially homogeneous models [15|, [IB]. Furthermore, since 
these integrability conditions are identically satisfied in the linearized theory (with 
FLRW background), Newtonian-like models have a linearization instability [|17], [T5|j , 
i.e. there are consistent linearized solutions which are not the limit of any consistent 
solutions in the full, nonlinear theory. 

The linearization instability found via an exact covariant analysis in [IT], [L5| is 
confirmed by an independent analysis based on a 1/c expansion of the field equations 
and Bianchi identities ||13|| . The gravito-magnetic tensor vanishes in the Newtonian, but 
not the post-Newtonian, limit. Forcing H a b = at the post-Newtonian level has the 
consequence that the non-local nature of the Newtonian tidal force cannot be recovered 
in the Newtonian limit, and reflects the fact that it is incorrect in general to assume 
H a iy = in general relativity. 

Thus the Newtonian-like silent models (H a b = 0) have a very narrow and limited 



applicability in cosmology and especially in gravitational instability [18], [15 
Realistic collapse scenarios or realistic inhomogeneous models of the late universe require 
a gravito-magnetic field EL a b- Spacetimes with H a b ^ do not have a Newtonian 
counterpart H, and they are consistent in the generic case, i.e. when there is also 



a gravito-electric field E a b and no external conditions are enforced on H a b or E a b [19 



The most 'extreme' non-Newtonian models (in terms of 'distance' from Newtonian 



§ A weaker covariant condition for the absence of gravitational radiation is that the spatial curls 
(defined below) of E a b and H a b must vanish 0, §, §] • 
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theory) are those with purely gravito-magnetic field, i.e. with E ab = 7^ H ab . (Note that 
E ab = and H ab = together imply that the spacetime is FLRW.) Consequently we 
will use the term 'anti-Newtonian' for irrotational dust spacetimes with E ab = 7^ H ab . 
Purely gravito-magnetic spacetimes (not restricted to irrotational dust) appear to have 



first been discussed in [20], where it was shown that either the shear or the vorticity must 
be non-vanishing. In this paper we consider only such models containing irrotational 
dust. Since E ab = also implies no gravitational radiation 0, the anti- Newtonian 
models are also 'silent'. 

Here we show that irrotational dust spacetimes with E ab = are subject to 
integrability conditions that are even more restrictive than in the H ab = case. The 
integrability conditions once again form in general non-terminating chains that lead to 
inconsistencies. There may be spatially homogeneous models of this type which satisfy 
the integrability conditions, but we have been unable to find examples. However, a 
further result indicates that there are unlikely to be any consistent exact solutions. 
This further result is that the only linearized irrotational dust solutions with E ab = 
are exactly FLRW (which have H ab = 0) - i.e. there are no linearized anti-Newtonian 
irrotational dust models. 

The question arises as to whether there are any known purely gravito-magnetic 
solutions (Eab = / H ab ) at all. In [^TJ, it is conjectured that there are no such 
non-flat vacuum solutions. An exact non- vacuum solution is given in [22], apparently 
the first purely gravito-magnetic solution. However, as pointed out in P2L this 
solution has unphysical Segre type. We show in an appendix that it is kinematically 
a magnetic counterpart of the Godel solution (which is purely gravito-electric 0), but 
that dynamically it has an unphysical source, with negative energy density. 

In section 2 we discuss briefly the covariant propagation and constraint equations, 
and the covariant approach that has been developed to analyze consistency We 
show in section 3 how the silent condition H ab = implies a primary covariant 
integrability condition, which is the basis for the linearization instability, as well as for 
the indefinite chain of conditions in the nonlinear case. This discussion forms a prelude 
to the analysis in section 4 of the E ab = case, where there are two primary integrability 
conditions, each of which produces in general an indefinite chain of further conditions 
after differentiation. In the linearized theory, the primary conditions themselves lead 
to the vanishing of anisotropy and inhomogeneity, i.e. to the FLRW case, so that there 
are no linearized anti-Newtonian models. 



We use the notation and conventions of 19, m. Units are such that 87rG = 1 



a,b, ■ ■ ■ are spacetime indices; (square) round brackets enclosing indices denote (anti-) 
symmetrization, while angled brackets denote the spatially projected, symmetric and 
tracefree part. 
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2. Covariant dynamical equations 

A covariant approach to the propagation and constraint equations and their consistency 
has been developed in [I9|. There it was shown that in generic irrotational dust 



spacetimes, i.e. without restrictions on E ab and H ab} the constraints are consistent 
and evolve consistently with the propagation equations. (See also |23|, El].) However, 
as discussed above and shown below, consistency breaks down in general when the 
Newtonian condition H ab = or the anti-Newtonian condition E ab = / H ab are 
imposed. 

The dust four- velocity u a (with u a u a = —1) provides a unique covariant 1 + 3 
splitting, as fully discussed in || . Here we follow the streamlined version of the formalism 
developed in 0, |J . The Weyl tensor splits covariantly into gravito-electric and gravito- 
magnetic partsf 

E a b — C ac bdU C U d , H ab = \ £ acdp Cd be 1 ^ > 

where the spatial permutation tensor is e a bc = Vabcd,u d , and r\ abc d is the spacetime 
permutation tensor. The tensor h ab = g ab + u a u b , where g ab is the metric tensor, 
projects orthogonal to u a . Then the projected tracefree symmetric part of a rank-2 
tensor is S^) = h a c h b d S^ c d) — \S c dh cd h ab . The fluid kinematics and dynamics, and the 
locally free field are described by the covariant scalars p (energy density) and (volume 
expansion rate), and rank-2 tensors a ab = a^ ab ) (shear), E ab = and H ab = H( ab y 

The pressure p and the vorticity u a are assumed to vanish in the irrotational dust case 
we consider here. 

The covariant derivative splits into a covariant time derivative S a ... = u b V b S a ... , 
and a covariant spatial derivative D iS&... = h c a h d ■ ■ ■ V c Sd--- ■ Then the latter leads to a 



covariant spatial divergence and curl [|I9| : 

div V = B a V a , curl V a = e abc B b V c , 
(div S) a = B b S ab , curl S ab = e cd (aD c S b) d . 

Important identities obeyed by these derivatives are collected in Appendix A. 

Covariant splitting of the Bianchi identities and the Ricci identity for u a , where 
the field equations are taken as an algebraic definition of the Ricci tensor, leads to the 
propagation equations 

p + Qp =0, (1) 

e + |0 2 = -\p-a ab a a \ (2) 

Oab + |0O"afe + 0"c(aO"b) C = ~ E ab , (3) 

E ab + QE ab - 3a c{a E b) c = curl H ab - \pa ab , (4) 



f See |5 



for a thorough discussion of gravito-electromagnetism. 
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H ab + 0H ab - 3a c{a H b) c = - curl E ab , (5) 

and the constraint equations 

C\ =DV a6 -|D a 6 = 0, (6) 

C 2 ab = curl a ab - H ab = , (7) 

C\ =D b E ab -lD aP -e abc a b d H cd = 0, (8) 

C 4 a ^D^ + e^a^^O. (9) 

The propagation equations (l)-(5) determine the covariant variables uniquely once 
initial data is specified on a surface S(t Q ), defined by t = to, where t is comoving 
proper time. However the constraint equations (6)-(9) place restrictions on the initial 
data, and must be satisfied on S(t) for all t. Since we have imposed the conditions p = 
(dust) and u a = (irrotational),f there is no a priori guarantee that the constraints will 
not lead to inconsistencies. However, consistency has been shown to hold in the generic 



case [[L9| . Lengthy tensor calculations lead to the following evolution equations of the 
constraints C A along u a : 

C\ = - QC\ + 2e a bc a b d C 2 cd - C s a , (10) 

C 2 a b = — QC 2 ab — ^ {a&fycC 1 d > (H) 

— — Ift/^S _|_ 1 bsi3 _ l n /7l 
^ O 3 WL -' a ' 2 a " 2 r a 

+ \E a h C\ - e a bc E b d C 2 cd + |curlC 4 a , (12) 

— — ift/ 1 ^ _|_ In- b r^ 

^ a 3 a ~r 2 " a ^ f) 

+ lH a b C\ - e bc H b d C 2 cd - icurlC 3 a . (13) 

It follows that if C A (t ) = 0, then C A = for all time t. Thus, if no further conditions 
(beyond p = and u a = 0) are imposed, the constraint equations are preserved 
under evolution in the generic case. We also require that the initial constraints are 
consistent, i.e. that C A (t ) = is not over-determined. We see this as follows [flQfl . 
If we freely specify, for example, a ab (to) and D o p(t ), then C 1 determines D o O(t ), C 2 
determines H ab (t ), and C 3 determines D 6 i? a f ) (to)- It could appear that C 4 then imposes 



a consistency condition, but this is not the case, since it can be shown that [19 



C 4 a = |curlC 1 -D 6 C 2 a6 . (14) 

This means that the constraint equation C 4 (t ) = is identically satisfied by virtue of 
C 1 (to) = = C 2 {to). Thus the constraint equations are consistent with each other, and 
they evolve consistently, in the generic case. 

f as well as the implicitly imposed conditions ir ab — (no anisotropic stress) and q a = (no energy 
flux) 
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Note that we can also interpret equation (]TJ|) as the statement that no new vector 
constraint arises from the divergence of the tensor constraint C 2 . In other words, the 
tensor constraint is essentially 'transverse traceless' in content. 

When additional covariant conditions are imposed, then this generally valid 
consistency may be disturbed and lead to integrability conditions. Consider the 
additional conditions that the gravito-electromagnetic tensors are divergence- free: 

B b H ab = 0^H ab and D b E ab = ^ E ab . 

These are covariant necessary conditions for gravitational radiation in linearized theory 
0, taken over into the nonlinear regime. Bianchi dust spacetimes were shown to include 
spatially homogeneous examples of such models in [IT5, and inhomogeneous G2 
solutions have also been found [p6 j. However, the Newtonian condition H ab = and 
the anti-Newtonian condition E ab = 7^ H ab both affect the propagation equations, 
converting one of them into a new constraint. This feature, which does not arise in the 
cases divif = or div E = 0, leads to complicated integrability conditions, as described 
in the following sections. 



3. Newtonian-like models (H ab = 0) 

When H ab = 0, the constraints C A are modified, but are still consistent, since equations 
(10)-(|14]) still hold. But there is now an additional constraint 

C 5 ab = curl E ab = , 

arising from the generic propagation equation (5), that must be satisfied, together with 
its evolution along u a . 

First we consider whether the divergence of the new tensor constraint C leads to 
an additional vector constraint. Using the identities (Al) and (A6), and the constraint 



equations (8) and (9), we find |T7| 

D b C\ b = |curlC 3 a - |6C 4 a - a a b C 4 b ■ (15) 

Thus no consistency condition arises from the divergence, and divC is determined by 
C A } where A = 1,2, 3, by virtue of (0). 

Now consider the evolution of C 5 . Using the identities (A5) and (A7), and the 



propagation equations (4) and (7), we get |T7 

' ab 



r 5 

ab 



4a/75 3_cd 771 n\ 

— ttUL ah — o£ („rJh\r^ t 



(a^b)cL 

\pC 2 a b ~ §£ Cd (aVfycC 3 d + \7~tab , 



(16) 



where 



K 



ab 



r -cd{ 



a {D e [E b) c a d e ] + 2D C [a b)e E de 
+ a b) c D e E d e + \a%\D e E h) d } . 



(17) 
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It follows that a necessary condition for consistent evolution of the constraints in 
Newtonian-like ('silent') universes is the covariant condition 

n ab = o. (is) 

This is the primary integrability condition for Newtonian-like models. Its repeated 



covariant time derivatives must also be satisfied. In |T5| , these derivatives, up to fourth 
order, are evaluated in a shear eigenframe, following the tetrad methods developed in 
27J, leading to a set of non-trivial conditions. In general, further derivatives produce 



further independent conditions, forming an indefinite chain of integrability conditions. 
The conditions are identically satisfied in Szekeres and spatially homogeneous silent 
models (including FLRW models, where E ab = 0), but not in general. Thus 
relativistic Newtonian-like models are in general inconsistent. The same result was 
found independently by different methods in ||16|| , and is supported by the further 
independent results of fl3f . 

Furthermore, equation (17) shows that the condition (|i~8l) and its time evolution are 
identically satisfied in the case of covariant linearization (see [^]) around an FLRW 
background characterized by 

D a p = D a 6 = and a ah = E ab = H ab = . 

Thus relativistic Newtonian-like models are subject to a linearization instability , in the 
sense that there exist consistent solutions of the linearized theory of this type which 
are not the limit of any consistent solution of the exact nonlinear theory. As pointed 



out in 1151 ) these results together cast serious doubt on the validity and usefulness of 
pursuing exact 'silent' solutions as realistic models of the late universe or of gravitational 
instability. Realistic general relativistic models involve a gravito-magnetic field, which 



is confirmed by the independent approach of [13]. However, as shown in the following 



section, a purely gravito-magnetic field leads to even more severe restrictions. 
4. Anti-Newtonian models (E ab = ^ H ab ) 

One of the nice features of Newtonian-like models, which facilitated extensive 
investigation of their dynamics, is that the condition = has the effect of decoupling 
the curls from the propagation equations. The propagation equations reduce to a 
coupled system of ordinary differential evolution equations, i.e. to equations (l)-(4) 
(with curl if a ft = in the last equation). Dynamical analysis of silent models based only 
on these equations involves the implicit assumption that the constraint equations are 



automatically satisfied. As shown in []T5j, [16] and outlined in the previous section, this 
assumption is incorrect. 

Interestingly, an analogous situation arises in the anti- Newtonian case E a b = ^ H a b 
(although, to our knowledge, these models have not previously been investigated). Once 
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again, as for silent Newtonian-like models, the curls decouple from the propagation 
equations, which again reduce to a coupled system of ordinary differential evolution 
equations. This feature reflects the fact that anti- Newtonian models are also ' silent '.f 

The anti-Newtonian propagation equations are given by (l)-(3) and (5), with 
E a b — in (3) and curl-E^ = in (5). The first 3 form a closed system determining 
the evolution of p, 6, and <j a b, with the last then determining the propagation of H ab 
without affecting the evolution of the other quantities. Thus the matter propagation 
is completely decoupled from the Weyl tensor. f Equation (4) with the left hand side 
vanishing is no longer a propagation equation - it has become a new constraint, as in the 
Newtonian-like models. One has to investigate the consistency of the new constraint, 
and this is done below. 

However, there is an added complication in the anti-Newtonian case, not present in 
the Newtonian-like models. The propagation equations are not independent, in the sense 
that propagation equation (5) for H ab must be consistent with the curl of propagation 
equation (3) for a ab , by virtue of the constraint equation (7). Using identity (A7), we 
can rewrite (5) as 

curl a ab + |9curl a ab - a e c e cd (aD e a b ) d = , 



and we find that the curl of (3) becomes 



curld ab + |6curlcr a6 + e cd ( a a b) d D e a ce + e cd{a D c 



<Jb)t 







where we also used identity (A5) and the constraint equation (6). The difference between 
these equations leads to the condition 



'cd(c 



a {D c [a de a b)e ] + D e [a b) d a 



0. 



This condition turns out to be an identity (satisfied by any tracefree symmetric tensor), 
which is given in [19|. Thus we can ignore the propagation equation (5) for H ab , since it 



follows from the shear propagation equation (3) and the constraint equation (6), using 
covariant identities. 

The coupled system of ordinary differential evolution equations to be satisfied is 
then 

P = -©P, (19) 

© = -|© 2 -|p-^ ab , (20) 

&ab = - f @cr a £> - a c ( a a b ) c . (21) 

f Note that within the class of perfect fluid spacetimes, the covariant propagation equations (see 
reduce to ordinary differential evolution equations under the more general conditions that ii a = and 
curl E ab = = carl H ab . 

f There is no Weyl tensor source term in the geodesic deviation equation, cf. Q. 
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These determine the evolution of the matter variables, which then determine the 
evolution of the gravito-magnetic field through 

H a b = —®H ab + 3cr c ( a H b ) c . (22) 

The constraint equation (9) shows that anti-Newtonian models have div H = 0. This 
is identically satisfied by virtue of equation ([14]), which continues to hold when E ab = 0. 
An additional constraint C 5 arises from the gravito-electric propagation equation (4). 
The system of constraint equations is then 

C\ = D b a ab -p a Q = 0, (23) 

C 2 ab = curiae -H ab = 0, (24) 

C\ = - ±D a p - s abc a b d H cd = , (25) 

C 4 a =D b H ab = 0, (26) 

C 5 ab ee curl H ab - \pa ab = . (27) 

We can eliminate H ab via the constraint equation (24), which expresses the fact that 
the shear is a covariant gravito-magnetic potential. Constraint equation (25) becomes 

e abc a b d cm\a cd = -§D a p. (28) 

Second-order derivatives then arise in constraint equation (27), which may be rewritten 

as 

D 2 a a6 = (±p - ±0 2 + a cd a cd ) a ab - Qa c {a a b)c + f D (a DV 6)c , (29) 

after using identity (A8) for the curl of the curl of a tensor, where D 2 = D a D a . 
The constraint equation (|29"D is a nonlinear generalization of the covariant Helmholtz 
equation. It may also be deduced as a special case of the nonlinear wave equation for 
the shear that is derived in fll9"| . 

The propagation equations (19)-(21) provide a unique solution for {p, 9, a a b}, given 
the values of these quantities on an initial surface S(t ). As in the Newtonian-like case, 
the problem is that the initial data is subject to a system of constraint equations [(23)- 
(27), or equivalently (23), ( p8[ ) and (0)], that is in general over-determined. Before 
showing this in detail, we can see intuitively how it arises, since the solution <J ab (to) of 
equation (|29|) , even though it allows for arbitrary tensors of integration, will in general 
not satisfy the constraint equation (f2~8|). 

We consider firstly whether the divergence of the new tensor constraint C 5 leads to 
an additional vector constraint. By the identity (A6) and the constraint equations (23), 
(25) and (9), we find 

D b C 5 afe = -yC\ + |0C 3 a + |curl C 4 a + \j a , (30) 

where 

J a = 0D Q p - 3pD a - WD bP . (31) 
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Thus by equation ([!(]) it follows that a necessary condition in anti-Newtonian universes 
for consistency of the constraints on an initial surface is the covariant condition 



pD a 9 = ±9D a p - la a b D bP . 



(32) 



This can be interpreted as an algebraic relation between the spatial gradients of p and G, 
or we can use the constraint equations (23)-(25) to rewrite it as an algebraic condition 
on the div and curl of the shear: 



9a ab e bcd a c e cur\ a de - 6Qe abc a b d cmlcr 



cd 







Equation (|32| ) is a primary integrability condition, whose successive derivatives must 
also be satisfied. There is at least one special situation where this condition is identically 
satisfied. If D a p = 0, as for example in spatially homogeneous models, then by the energy 
conservation equation (19) and the identity (A2), we find D a = 0. It then follows that 
J a is identically zero. However, note that the new constraint equation (27) itself is not 
necessarily identically satisfied when ~D a p = - only its divergence vanishes identically, 
as seen from equation (|30|) . Below we will find a further primary integrability condition 
from (27) which is not identically satisfied when D a p = 0. We have been unable to 
find spatially homogeneous solutions that satisfy equation (27) when H ab ^ 0, and 
the existence of such solutions remains an open question. (Clearly FLRW solutions, 
with H ab = = o~ a b, satisfy equation (27) identically.) Indeed, the only solution with 
E a b = 7^ H ab known to us is not an irrotational dust solution (see Appendix B). 

In general, and especially in the more physically interesting cases, D a p is nonzero 
and the condition ( |32|) is not trivial. The evolution of integrability condition (^) along 
u a produces a further integrability condition. Using identity (A2) to commute time and 
space derivatives, propagation equations (19)-(21) to eliminate time derivative terms, 
and condition (K3) to eliminate D a 0, we get 



7~2\ b 









-§p + U° 2 )A D «P - i®°abV b p - \(a 2 ) {ab) D b p, (33) 



where (a 2 ) ab = a a c a bc is the contracted tensor product. It is apparent that in general, 



with D a p 7^ 0, condition (|33D is not automatically satisfied if (|32j) is, i.e. the derived 
condition (|33|) is not an automatic consequence of the primary condition (p2|). A further 
time derivative gives 



pD a 



§6(aV + W)A d *p 

- i(a 2 ) c c a ab B b P - |6(a 2 ) (afc) DV - ^) {ab) D b p , (34) 



where we used the conditions ( |3"2]) and 
an integrability condition of the form 



). Clearly the iV-th time derivative leads to 



'-N+l\ b 



A [N+1) D a p + A {N) tT ab D b p + ■■■ + A (0) (a N+1 ) {ab) D b p 
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where A^m) involves in general p, and (a 1 )^, I — 0, 1, • • • , M. The algebraic invariants 
of the shear are {<J 2 ) a a and (a 3 ) a a (the maximal three further independent components 
of the shear correspond to the rotational freedom in the choice of frame). The scalars 
(v N )a a are not all independent. 

Thus there is an indefinite chain of derived integrability conditions on S(to), all of 
which must be satisfied. At each level, the condition does not follow automatically from 
lower-level conditions. Since each such equation involves only the initial data {p, 6, o~ ab } 
on S(t ), it is clear that in general the chain of conditions is over-determined and will lead 
to inconsistencies, i.e. relativistic anti-Newtonian universes are in general inconsistent. 
We conjecture that the new constraint equation (27) and the integrability conditions 



(p2|), (p3|) , . . . that follow from it are only consistent if the shear, and hence the gravito- 
magnetic field, vanishes. However, we have been unable to prove this conjecture given 
the complicated nature of the chain of integrability conditions. 

This is similar to the situation in silent Newtonian-like models, but in contrast to 
the case of shear-free rotating dust models. Using tetrad methods as opposed to the 
covariant approach adopted here, a condition similar to (p2|) - i.e. an algebraic relation 
between spatial gradients - was found for shear- free dust in p7| , and its successive time 



derivatives led to the conclusion that Qu a = 0. Recently, this result was regained via 
covariant methods in [^UJ (and then generalized in |n| from dust to vanishing fluid 
four-acceleration). The crucial simplifying factor in the shear-free case is the central 
role of the naturally defined vector u a , as opposed to the irrotational case, where there 
is no natural algebraically defined vector, and instead the central role is played by the 
tensor o ab . Furthermore, the propagation equation for oj a is linear in cj a , in contrast to 
the shear propagation equation (21), which is nonlinear in a ab , and which gives rise to 
the proliferation of terms involving the gradient of the trace of tensor products. 

The conclusion arising from the spatial divergence of the anti-Newtonian constraint 
equation (27) is that in general the models are not consistent. The conjecture is that 
there are no consistent anti-Newtonian models. Both of these are reinforced by the 
existence of a further chain of integrability conditions. This arises from the time 
evolution of (27), which must also be satisfied. Using the identities (A5) and (A7), 
the propagation equations (19) and (21), and the constraint equations (23) and (24), 
we find thatf 



C 5 ab — — |@C 5 a fe — §£ Cd (aHtycC 1 d + S ab , (35) 



where 



o~ c ( a Hb) c 



Sab = \p®o- ab + \pa c{a a b) c + 3H c{a H b) c + 3curl 

V e H b) d . (36) 



f Instead of constraint equation (24), we can also use the gravito-magnetic propagation equation (22). 
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It follows that a necessary condition for consistent evolution of the constraints in anti- 
Newtonian universes is the covariant condition 

£ab = 0. (37) 



Note that, in contrast to the previous integrability condition (p2|), this condition involves 
second derivatives of the shear, given that H ab = curler^. 

Clearly this condition is identically satisfied in the degenerate case of a a b = 0, but, 
in line with the conjecture stated above, it is unlikely to be satisfied for any shearing 
solutions because of the chain of further conditions that are implied. Specifically, the 
evolution of this integrability condition must also be satisfied. As before, the propagation 
equations can be used to eliminate time derivatives and arrive at a chain of derived 
integrability conditions intrinsic to S(to), which places inconsistent restrictions on the 
initial data {p, O, a a b} in addition to those remarked on above. Since inconsistency of 
the anti-Newtonian models in general is already implied by these previous restrictions, 
we will not give the very complicated condition arising from S a b = which, together 
with its time derivatives, strongly reinforces the conclusion arrived at already. 

The key role of the condition (^) emerges in the linearized case. Our conjecture 
that there are in fact no consistent anti-Newtonian models is further reinforced by an 
examination of the linearized form of the integrability conditions. Linearization about an 
FLRW universe (see ^8[) of (|37D shows that, in contrast to the Newtonian-like case, the 



linearized integrability condition is non-trivial, so that not all linearized anti-Newtonian 
solutions are consistent. In fact none are consistent, since the linearized integrability 
conditions are satisfied only when H a b = 0. Since already E a b = 0, this is the FLRW 
case, which is Newtonian-like. This result can be derived as follows. The linearization 
about an FLRW background of (|37D produces 

ea afe = 0, (38) 

implying either O = or a a b = 0. The linearized form of the other primary integrability 
condition (|32T) is 



3pD a 6 - 0D a p = , 

which is automatically satisfied if G = 0, and also holds if a a b = 0, since in that 
case the spacetime is FLRW. However FLRW spacetimes are not anti-Newtonian, and 
= implies via the linearization of propagation equation (20) that p = 0, and we 
have ruled out this vacuum case by our definition of dust universes. (Furthermore, the 
exact nonlinear form of propagation equation (20) shows that = 0^ a a b leads to the 
unphysical condition p < 0.) 

Thus there are no linearized anti-Newtonian universes. The implication of this result 
is that it is difficult to see how any consistent exact anti-Newtonian solution can exist. 
Such a solution would need to have the property that it cannot be linearized about an 
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FLRW solution - the solution could not admit small gravito-magnetic field or shear. If 
such solutions existed, the FLRW solution would have to be an isolated point in the 
space of all irrotational dust solutions with E ao = 0. 



5. Concluding remarks 

In summary, the overall conclusion following from the covariant analysis in sections 3 
and 4 is that irrotational dust universes in general relativity with realistic inhomogeneity 
must have both gravito-electric and gravito-magnetic fields. The Newtonian-like case 
H ab = is too restrictive, supporting the argument that there is not a straightforward 
relationship between general relativistic and Newtonian universes, while the anti- 
Newtonian case E a b = ^ H ao is even more restrictive. Newtonian-like universes 
(i.e. those with = 0) are in general inconsistent, and subject to a linearization 
instability. We showed that anti-Newtonian models are also in general inconsistent 
by virtue of the indefinite chain of integrability conditions fl3"2|), (|33]), . . . , which arise 
from the spatial divergence of the new tensor constraint, and which over-restrict the 
shear. We conjectured that only the degenerate shear-free case, i.e. the FLRW models, 
satisfy these integrability conditions. This conjecture is strongly reinforced by the second 
integrability condition (|37|), which arises from the time evolution of the new constraint, 
and whose linearized form implies Q<J a b — 0, leading to the non-existence of any anti- 
Newtonian solutions linearized about FLRW models. 

These results extend our previous work on consistency within the class of irrotational 
dust spacetimes |0|, [H| [23) , using a streamlined version |j] of the covariant 1 + 3 



formalism ||, and a systematic covariant approach to analysing consistency [19|. This 



approach parallels the tetrad methods developed in |27|] for analysing the consistency of 



shear-free dust spacetimes, and is well adapted for investigating more general classes of 
spacetime. 
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Appendix A. Covariant Identities 

For convenience we collect here the necessary identities from [TP 



curl D a f 
(Da/)' 

(B a v b y 
(n b s ab )' 

curl (fS ab 
D b curl S, 



ab 



(curl S^)' 
curl curl S ab 



0, 

Dj-lQD a f-a a b D b f, 

D a V b - |0D a H - a a c D c V b + H a d e dbc V c 

D b S ab — ^QD b S ab — a bc D c S ab + e abc 

y abc D b (D d S cd ) + e abc S b d (l&a cd - E 
— cr a b£ a ce D d , 



■H b d S cd 



red 



'ab 



curl S ab — -|0curl S, 

— Ce C £cd(aD e S b ) d + 

— £> 2 s ab + |D< a D c s&) 



— &e C £cd(aD e S b ) d + 3H c ( a S b )° , 



E b ) c — \QcF b ) 



'ab 



— ° C a<2b d Scd 



>b) 

S cd + cr c ( a S b ) a cd 



+ o ab a cd S, 



cd 



where S ab = S( ab ) and D 2 = D a D a is the covariant Laplacian. 
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Appendix B. A purely gravito-magnetic solution 



In P2| , an exact solution of Einstein's field equations with E a b = 7^ H a b is found, 
apparently the first such solution. The solution is given in the form (reversing the 
signature to conform with our convention) 

ds 2 = - x~ 2 du 2 + (l - x 2 ) dy 2 - 2du dy 

+ x~ 2 exp (-Ax 2 ) dx 2 + x~ 2 dv 2 , (B9) 



where we have labelled the coordinates in the sequence x a = (u,y,x,v). The four- 
velocity u a is proportional to l a + n a (there is a misprint in |2!| that turns the + into 
a — ), where 

I = 2~ 1/2 [-zrMw + (1 - x)dy] , n = 2~ 1/2 [-x^du - (1 + x)dy 

are Newman-Penrose null vectors (using the notation v = v a dx a ). The Segre type of 
the nonzero Ricci tensor is given in |22| as {1 1 z z}, but no further discussion is given 
of the properties of (B9). 
We find that 

u = —x du — xdy , 

which implies 

u A du = 2x~ 1 du A dx A dy ^ . 

It follows that the solution is rotating. Furthermore, it is apparent from (B9) that 
£ a = 5 a o is a timelike Killing vector, so that the solution is stationary. (It is not 
static, since £ A d£ 7^ 0, or, equivalently, since £ a is parallel to u a .) Since u a = x£ a , 
Killing's equation shows that the solution is non-expanding and non-shearing. Thus it 
is kinematically characterized by 

6 = , a ab = , u a + . (BIO) 

This solution therefore has similar kinematic characteristics to the Godel solution, and 
we can think of it as a gravito-magnetic counterpart of that solution, which is purely 
gravito-electric {ti a b = 0) because V&u; a = 0.f 

However, unlike the Godel solution, the magnetic solution (B9) has no physically 
significant interpretation by virtue of its pathological Segre type (see pifl , p. 72). For 



this Segre type of the Ricci tensor, the weak energy condition is violated, so that 
the energy density is negative. Not only does the solution fall outside of the class 
of irrotational dust models, but it also has an unphysical source. 



f Note that the Godel solution provides a counter-example to the notion that rotating matter always 
produces a gravito-magnetic field . 
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